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We study the resonant spin-flavor (RSF) conversion of supernova neutrinos, which is induced by 
the interaction between the nonzero neutrino magnetic moment and the supernova magnetic fields, 
and its dependence on presupernova models. As the presupernova models, we adopt the latest 
ones by Woosley, Heger, and Weaver, and, further, models with both solar and zero metallicity are 
investigated. Since the (1 — 2Y e ) profile of the new presupernova models, which is responsible for the 
RSF conversion, suddenly drops at the resonance region, the completely adiabatic RSF conversion 
is not realized, even if \i v Bq = (10 _12 /xb)(10 10 G), where Bo is the strength of the magnetic field at 
the surface of the iron core. In particular for the model with zero metallicity, the conversion is highly 
nonadiabatic in the high energy region, reflecting the (1 — 2Y e ) profile of the model. In calculating 
the flavor conversion, we find that the shock wave propagation, which changes density profiles 
drastically, is a much more severe problem than it is for the pure Mikheyev-Smirnov-Wolfenstein 
(MSW) conversion case. This is because the RSF effect occurs at a far deeper region than the MSW 
effect. To avoid the uncertainty concerning the shock propagation, we restrict our discussion to 
0.5 s after the core bounce (and for more conservative discussion, 0.25 s), during which the shock 
wave is not expected to affect the RSF region. We also evaluate the energy spectrum at the Super- 
Kamiokande detector for various models using the calculated conversion probabilities, and find that 
it is very difficult to obtain useful information on the supernova metallicities and magnetic fields or 
on the neutrino magnetic moment from the supernova neutrino observation. Future prospects are 
also discussed. 

PACS numbers: 95.85.Ry, 13.40.Em, 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

A core-collapse supernova explosion is one of the most 
spectacular events in astrophysics; 99% of its gravita- 
tional binding energy is released as neutrinos, and only 
1% as the kinetic energy of the explosion. Therefore, neu- 
trinos play an essential role in supernovae, and their de- 
tection by ground-based large water Cerenkov detectors, 
such as Super-Kamiokande (SK) and Sudbury Neutrino 
Observatory (SNO), would provide valuable information 
not only on supernova physics but also on the nature 
of neutrinos. What we can learn from the next galac- 
tic supernova has been considered in many articles (for 
a review, see Ref. Q). For example, we can constrain 
the properties of neutrino oscillations, such as the mixing 
angle between the first and third mass eigenstates (#13), 
and the mass hierarchy [normal (mi <C 771,3) or inverted 
(mi »m 3 )] Hi. 

In addition to the nonzero neutrino masses and mix- 
ing angles, the nonzero magnetic moment is of a diffcr- 
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ent nature for neutrinos beyond the standard model of 
particle physics, and has attracted a great deal of atten- 
tion from many theoretical and experimental physicists. 
If neutrinos have a nonzero magnetic moment, it leads 
to precession between left- and right-handed neutrinos 
in sufficiently strong magnetic fields [3-0] In general, 
nondiagonal elements of the magnetic moment matrix 
are possible and neutrinos can be changed into different 
flavors and chiralities Furthermore, with the ad- 

ditional effect of coherent forward scattering by matter, 
neutrinos can be resonantly converted into those with 
different chiralities 0, 0, 0] by a mechanism similar to 
the well-known Mikheyev-Smirnov-Wolfenstein (MSW) 
effect [HI O This resonant spin-flavor (RSF) con- 
version induced by the neutrino magnetic moment in 
strong magnetic fields was first introduced to solve the 
solar neutrino problem, and actually gav e the best fit so- 
lution before the KamLAND result [l^. However, the 
recent KamLAND experiment 0] has shown that the 
large mixing angle MSW solution is the most favorable 
one; the RSF mechanism is suppressed at the subdomi- 
nant level. From the KamLAND negative results for the 
solar antineutrino search, an upper bound on the neu- 
trino magnetic moment is obtained, /j^lx 10~ 12 /is, 
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where is the Bohr magneton |16(. This upper bound 
is comparable to the most stringent limit from the stellar 
cooling argument, fi u <(1-4)x10 _12 /xb |l7j . 

Although the RSF mechanism does not work at a 
dominant level in the Sun, it may occur efficiently in a 
denser environment with stronger magnetic field, which 
is actually expected in the case of core-collapse super- 
novae. The RSF conversion mechanism in supernovae 
has been investigated by many authors [el 1^. [Tot fl8t liH 
El El SIliliEii. Among them, Ando and Sato 
[26j| studied the RSF effect using a three-flavor formula- 
tion with the latest oscillation parameters, and pointed 
out that the combination of the MSW and RSF effects 
makes the crossing scheme very interesting to investigate. 
They showed that until 0.5 s after the core bounce the 
RSF-induced v e <-*■ v r transition occurs efficiently, when 
fJtu £ f(r 12 ^B(£ /5xf0 9 G) _1 , where B is the strength 
of the magnetic field at the surface of the iron core. 

The effective matter potential for the RSF conversion 
is given in a form proportional to the value of (1 — 2Y e ), 
where Y e is the electron number fraction per nucleon. 
Thus, the deviation of the value of Y e from 0.5 in the stel- 
lar envelope is quite important, and this value is strongly 
dependent on the isotopic composition. Since this devia- 
tion is determined by rarely existent nuclei, an accurate 
estimate of this deviation is quite difficult. Therefore, 
the astrophysical uncertainty in (f — 2Y e ) should be dis- 
cussed. This point was first investigated by Totani and 
Sato [23 . However, their treatment was based on a two- 
flavor formulation with the large uncertainties concerning 
the mixing parameters of those days. Therefore, we need 
definitive investigations using the three-flavor formula- 
tion with the latest presupernova models as well as the 
recently determined neutrino parameters. 

In this paper, we study the RSF conversion mecha- 
nism using the three-flavor formulation with the latest 
neutrino mixing parameters. In particular, we investi- 
gate the dependence on presupernova models; we use the 
latest 15-Mq model by Woosley et al. [27|. and compare 
the results with those obtained with the previous pro- 
genitor model by Woosley and Weaver |28|. which was 
also adopted in the calculations of Ando and Sato [26| . 
It is also expected that the value of (1 — 2Y e ) strongly 
depends on the stellar metallicity, and hence we use the 
15Mq model with two different metallicities, solar and 
zero metallicity, and the metallicity effect is investigated. 

Throughout this paper, we adopt the realistic neu- 
trino mixing parameters inferred from the recent experi- 
mental results: for the atmospheric neutrino parameters, 
Am 2 3 = 2.8 x 10~ 3 eV 2 , sin 2 26» 23 = 1.0, and for the solar 
neutrino parameters, Am 2 2 = 5.0 x 10~ 5 eV 2 ,tan 2 612 — 
0.42. As for the still uncertain parameter 6*13, we as- 
sume sin 2 2^i3 = 10 -6 , 1 and we also assume a normal 



mass hierarchy. For the elements of the neutrino mag- 
netic moment matrix ji^j , where i and j denote the flavor 
eigenstates of the neutrinos, i.e., e,/i, and r, we assume 
that all the values of fiij are near the current upper limit, 
or fiij = 10~ 12 [i B - 

The remainder of this paper is organized as follows. In 
Sec. [n] we give the formulation used in our calculation, 
which includes all three-flavor neutrinos and antineutri- 
nos, and from the formulation a level crossing diagram, 
which enables us to understand the conversion scheme 
intuitively, is introduced. In Sec. IIIII the presupernova 
model we adopt in our calculations is illustrated, and 
the results of numerical calculations are shown in Sec. 
IIVI Finally, detailed discussions of the model uncertain- 
ties and of whether we can obtain, at present or in the 
future, useful information on physical and astrophysical 
quantities are presented in Sec. 

II. FORMULATION AND LEVEL CROSSING 
SCHEME 

A. Interaction with matter and magnetic fields 

The interaction of the magnetic moment of neutrinos 
and magnetic fields is described by 

((^)fll^intl(^)i) =(-LijB ± , (1) 

where B± is the magnetic field transverse to the direc- 
tion of propagation, and {v)r and {v)l are the right - 
and left-handed neutrinos, respectively. If neutrinos are 
Dirac particles, right-handed neutrinos and left-handed 
antincutrinos are undetectable (sterile neutrinos), since 
they do not interact with matter. On the other hand, 
if neutrinos are Majorana particles, i/r's are identical to 
antiparticles of vl§ and interact with matter. In this 
paper, we assume that neutrinos are Majorana particles. 
The diagonal magnetic moments are forbidden for Ma- 
jorana neutrinos, and therefore only conversion between 
different flavors is possible, e.g., (v e )R (^,t)l- 

Coherent forward scattering with matter induces an 
effective potential for neutrinos, which is calculated using 
weak interaction theory. The effective potential due to 
scattering with electrons is given by 

V±± =±V2G F (±^ + 2sm 2 dw^ n e , (2) 

where n e is the electron number density, Gf is the Fermi 
coupling constant, and 9w is the Weinberg angle. The 
± sign in front refers to v (+) and v (— ) and that in 
the parentheses to v e (+) and !/„ )T (— ). The difference 
between e and fi, t neutrinos comes from the existence 
of charged-current interaction. The subscript ±± of V 



1 As already shown in Ref. |2(1 . the dependence of the detected 
event on the parameter 813 is very weak, and even if we assume 



a large value of 613, the results do not change. 
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refers to the first and the second ± sign. The ordinary 
MSW effect between v e and z/„ T is caused by the poten- 



tial difference V e — V^ T = V++ — V-\ = \/2Gpn e . To 



(U,T 

+ 



include the RSF effect, which causes conversion between 
neutrinos and antineutrinos, we should take into account 
the neutral-current scattering by nucleons: 



where p is the density, is the nucleon mass, and 
Y e = n e /(n e + n n ) is the number of electrons per baryon. 
[When we obtained Eq. (@J, we assumed charge neutral- 
ity n e = n p .} 



V = V2G F [ I - 2 sin 2 9 W 



(3) 



where n p , n n are the proton and neutron number density, 
respectively. For neutrinos we add +V to the potential 
and for antineutrinos —V. Therefore, the RSF conversion 
between 9 e and Vu,t obeys the potential difference 



Ay = v ? -v, 



(V-+ -v)- (V+ 
V2G P 



V) 



r F (f 

m N 



(4) 



B. Three-flavor formulation 

The time evolution of the flavor mixing, which includes 
all three-flavor neutrinos and antineutrinos, is described 
by the Schrbdinger equation 



dr 



H B±M 
-B X M H 



(5) 



where 







#o = ( Am( 2 | C/ T + ( V+-+V | , (7) 

V+-+V 

- V 

H ° = 7^r U I Am 12 | C/ T + ( - V I , (8) 

V—-V 

U e l U e2 U e3 \ / C12C13 S12C13 S13 

U = \ Uf+i U,j,3 = -S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 ) , (9) 

U T i U r2 U t3 I \ S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 



M = 




(10) 



and Cij — cos 6>y , Sjj = sin 6ij . [We assume the CP phase 
S = in Eq. J^J for simplicity] 

Resonant flavor conversion occurs when two diagonal 
elements of the matrix in Eq. JSJ have the same value. 
There are five such resonance points, which are for v e «-> 
Up (MSW-L), u e <-> v T (MSW-H), v e (RSF-L), 
v e <-* v T (RSF-H), and «-* v T conversions. The suffixes 
"-L" and "-H" attached to "MSW" and "RSF" indicate 
whether the density at the resonance points is lower or 
higher. Hereafter, we neglect the <-> v T conversion, 
since it is always nonadiabatic and including it further 
complicates the discussion. 



Figure ^ shows the level crossing diagram, which we 
introduced in our previous paper 26] , to understand the 
flavor conversions described by Eq. JSJ intuitively. The 
figure clearly includes not only the ordinary MSW reso- 
nances but also the RSF effects, and it is expected that 
the combined effect of MSW and RSF makes this scheme 
very interesting to investigate. For instance, when MSW- 
L and RSF-H are adiabatic and the others are nonadi- 
abatic (this case is actually expected if #13 is small and 
the magnetic field is sufficiently strong), conversions such 
as v e -> ia>,^ -> V\,v' r -> Ux,v e -> jaj,!^ -> i> 2 , and 
v' r — > D3 occur. We can easily predict this sort of con- 
version scheme from Fig. ^ when all the resonances are 
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FIG. 1: Schematic illustration of level crossings, where ^1,2,3 
and ^1,2,3 represent the mass eigenstates of neutrinos and an- 
tineutrinos in matter, respectively, and v'„ T and pf, T the mass 
eigenstates at production, which are superpositions of and 
v T or and v T . There are four resonance points, MSW- 
L, MSW-H, RSF-L, and RSF-H. Adiabatic conversion means 
that the neutrinos trace the solid curve at each resonance 
point (i.e., the mass eigenstate does not flip), while nonadia- 
batic conversion is shown by the dotted line. 



either completely adiabatic or completely nonadiabatic; 
for the intermediate cases we have no choice but to trust 
numerical calculations. 
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FIG. 2: Presupernova profiles (W02S) used in our calcula- 
tions taken from Ref. |27|. Upper panel: The density and 
Y e combination that is responsible for the RSF conversions 
[|p(l — 2F e )|, solid curve], and that for the MSW conversions 
(pYe, dashed curve). Two horizontal bands represent A12 
and A 13 (these definitions are given in the text); at the inter- 
sections between them and the profile curves, the RSF and 
MSW conversions occur. Lower panel: The mass fraction of 
the various elements. 



III. PRESUPERNOVA MODELS 
A. Density and Y e profiles 

In this subsection, we discuss the several presupernova 
models that we adopt in our transition calculations. In a 
previous paper 26| , we used only the presupernova model 
of Woosley and Weaver [2^], which is for 15M Q and so- 
lar metallicity (hereafter W95S, where "S" denotes so- 
lar metallicity). However, we should investigate the de- 
pendence on adopted presupernova models with various 
metallicities. This is because the RSF conversion is very 
sensitive to the deviation of Y e from 0.5 [see Eq. (0}], 
which strongly depends on the metallicities as well as on 
the weak interaction rates adopted in the simulation of 
stellar evolution. 

In this study, we adopt the latest presupernova mod- 
els by Woosley et al. [27] with both solar and zero 
metallicities: W02S and W02Z ("Z" denotes zero metal- 
licity). Figures [21 and show density profiles (up- 
per panel) or pY e , which is responsible for MSW, and 
\p(l — 2Y e )\ for RSF, and the composition of each el- 
ement (lower panel), for the W02S and W02Z models, 
respectively. In the upper panels of both figures, we also 
show A12 = mjq /S.m\ 2 cos 29 i2/2^/2G fE v and A13 = 




FIG. 3: The same as Fig. H but for the W02Z model. 
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rriN AtiJ] 3 cos 2Q\^I2^2GfE v as two horizontal bands 
(the bandwidth comes from the energy range 5-70 MeV). 
At intersections between A12, A13 and p(l— 2Y e ), pY e , the 
RSF and MSW conversions take place. 

For the W02S model, the region where RSF-H occurs 
is the silicon burning shell, in which (1 — 2Y e ) suddenly 
drops, in contrast with RSF-L, where p(l — 2Y e ) gradu- 
ally changes. On the other hand, for the W02Z model, 
(1 — 2Y e ) suddenly becomes exactly zero at the bound- 
ary between Si + O and O + Ne + Mg layers, and this 
tendency continues to the He layer because of the lack of 
heavy nuclei, which cause the deviation of Y e from 0.5. 
In consequence, both the RSF-H and RSF-L conversions 
are expected to be highly nonadiabatic, as discussed in 
more detail in Sec. IIVI Note that the pY e profiles of the 
two different models agree well with each other, which 
clearly indicates that the ordinary MSW conversions are 
not noticeably affected by the value of Y e . 




10-= r -j 
10-' r -j 

10-= E r 
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FIG. 4: The \p(l — 2Y e )\ profiles, which are responsible for 
the RSF conversions, for the W02S (solid curve) and W95S 
(dashed curve) models. 

In Fig. 01 we compare the \p(l — 2Y e )\ profile of the 
W02S model with that of the W95S model, which was 
adopted in previous publications including Ref. |26( . The 
difference comes from the included weak interaction rates 
for nuclei. In the latest model W02S, a recent shell model 
is included in the calculations and results in substantial 
revisions to the older data set in the W95S model (see 
Ref. [2^1 for a detailed discussion). As a result, the value 
of (I — 2Y e ) decreases by a few orders of magnitude be- 
tween the He and Si + O shells. Particularly with the 
new model W02S, there is a sudden drop of p(l — 2Y e ) 
when the RSF-H conversion occurs, in contrast with the 
gradual decrease in the case of the W95S model. Since 
the adiabaticity of the resonance is reciprocal to the gra- 



dient of the logarithmic value of p(l — 2Y e ), the RSF 
conversion is expected to be less efficient than that with 
W95S given in Ref. ||. 

Although we use static progenitor models in calcu- 
lating the flavor transition, in fact the density profile 
changes drastically during a neutrino burst (~ 10 s) ow- 
ing to shock wave propagation, and we should use the 
time-dependent profiles 30] . Unfortunately, however, su- 
pernova explosion mechanisms are still controversial, and 
there is no reliable model that precisely describes the 
time-dependent density and Y e profiles. Further, there 
is also a large uncertainty concerning the magnetic field 
structure which is affected by the shock wave propaga- 
tion. Therefore, from this point on, we confine our dis- 
cussion to 0.5 s after core bounce, since in that case using 
the static presupernova and magnetic field models is con- 
sidered to be a good approximation [2(|. This is based 
on the numerical calculation by Takahashi et al. [3l| , the 
only authors having succeeded in shock propagation to 
the outer envelope. Although we cannot trust the details 
of their result without any doubt, the choice of the time 
scale during which the shock effect can be neglected is 
expected to be reasonable. In addition, for more conser- 
vative discussions, we also give calculations during first 
0.25 s. 



B. Magnetic fields 

We assume that the global structure of the magnetic 
field is a dipole moment and the field strength is nor- 
malized at the surface of the iron core with the values 
of 10 8 and 10 10 G. The reason for this normalization is 
as follows. The magnetic fields should be normalized by 
fields that are static and exist before the core collapse, 
because those of a nascent neutron star can hardly af- 
fect the far outer region, where the RSF conversions take 
place, within the short time scale of a neutrino burst. 
As discussed in the previous subsection, since the shock 
wave does not affect the resonance region at < 0.5 s af- 
ter bounce, it is also expected that the magnetic field 
structure and strength at the resonance points are not 
seriously changed at that time. The strength of such 
magnetic fields above the surface of the iron core may be 
inferred from observations of the surface of white dwarfs, 
since both are sustained against gravitational collapse by 
the degenerate pressure of electrons. Observations of the 
magnetic fields in white dwarfs show that the strength 
spreads in a wide range of 10 7 -10 9 G Considering 
the possibility of the decay of magnetic fields in white 
dwarfs, it is not unnatural to consider magnetic fields up 
to 10 10 G at the surface of the iron core. Then, in Eq. 
©, B± = B (r /r) 3 sin 9, where B is the strength of 
the magnetic field at the equator on the iron core surface, 
ro the radius of the iron core, and O the angle between 
the pole of the magnetic dipole and the direction of neu- 
trino propagation. Hereafter, we assume sin 6 = 1. 
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IV. RESULTS 



A. Conversion probabilities 



sin 2 2(9 13 =10- 6 ; B =10'°G 




o.oi o.i 

Radius [R ] 



magnetic field Bq = 10 10 G, the — > v e transition is 
almost completely adiabatic over the entire energy range. 

For all these three models, it appears that the RSF-L 
conversion does not play any role. This is because for 
the W02S and W95S models RSF-L occurs farther out 
than RSF-H, where the magnetic field strength is not 
large enough to induce adiabatic conversions. On the 
other hand, for the W02Z model, the RSF-L occurs in 
the same place as RSF-H; the magnetic field strength 
is also the same for these two resonances. In this case, 
however, the very abrupt drop of the p(l — 2Y e ) profile 
at RSF-L strongly suppresses efficient flavor conversions. 

At the end of this subsection, we focus on the decrease 
of the conversion probability at > 0.1 Rq which can be 
seen in Fig. [SJ This does not indicate any resonances 
but is merely an effect of flavor mixings. At such a deep 
region in the supernova envelope as < O.IRq, the P e 's 
propagate like the mass eigenstates owing to the large 
matter potential; the v e 's do not mix with the other fla- 
vor antineutrinos in this region. As they propagate to 
a radius that is larger than ~ O.li?©, the matter poten- 
tial becomes smaller; at this stage the D e 's are not the 
mass eigenstates at all, and mixing with the other flavor 
antineutrinos occurs. This effect induces the decrease at 
> 0.1i? Q seen in Fig. |S| 



FIG. 5: Conversion probability P{u^ — > v e ) as a function of 
radius for Bo = 10 10 G. The probabilities calulated with the 
W02S (a), W02Z (b), and W95S (c) models are plotted. 



We calculated Eq. (JHJ numerically with the adopted 
models given in Sec. II I II and obtained the conversion 
probabilities for each flavor. Among them, we show in 
Fig. [5] those of the Vy, — > v e transition for various presu- 
pernova models in the case of Bq — 10 10 G. (This conver- 
sion channel is essential in order to discuss the efficiency 
of the RSF effects.) 

For the W02S model, Fig. shows that the conver- 
sion between f MjT and v e occurs at a radius independent 
of energy, and becomes adiabatic as the energy increases; 
this character reflects the p(l — 2Y e ) profile given in Fig. 
121 which shows a sudden drop in the RSF-H region. 

The behavior changes dramatically when the metal- 
licity is zero, or for the W02Z model. In this case, as 
shown in Fig. E{b), the transition is completely adia- 
batic at low energy; however, once the energy is increased 
beyond some critical value, the conversion abruptly be- 
comes almost completely nonadiabatic. This tendency 
also reflects the characteristic profile in Fig. [3] 

The flavor transition occurs most efficiently in the case 
of the W95S model as shown in Fig. [SJc), which indicates 
the most moderate profile at each resonance point (see 
the dashed curve in Fig. With a sufficiently strong 



2 Actually, the results for the W95S model are taken from previous 



B. Energy spectrum at the Super-Kamiokande 
detector 

With the conversion probabilities given in the previ- 
ous subsection and the original neutrino spectrum emit- 
ted from the supernova core, we can calculate the flux 
of each flavor neutrino. From this point on, we assume 
that the distance to the supernova is 10 kpc. As the 
original neutrino spectrum, we adopt the results by two 
groups: the Lawrence Livermore group |33| and Thomp- 
son et al. |23. The Livermore spectrum [3i| resulted 
from a calculation with 20M Q progenitor models, and 
we label it LL20. Recently, that calculation has been 
criticized, since it lacks relevant neutrino processes such 
as neutrino bremsstrahlung and neutrino-nucleon scat- 
tering with nucleon recoils, which were not recognized to 
be important at date of calculation. However, since there 
are no other successful simulations of a supernova explo- 
sion, it is premature to conclude that their result is no 
longer reliable and we adop t their results. On the other 
hand, Thompson et al. [34J calculated for three different 
mass progenitors 11M©, 15M®, and 20M Q , and we label 
the models TBP11, TBP15, and TBP20, respectively. 
Although they did not succeed in simulating the explo- 
sion and their data end at 0.25 s after the core bounce, 



calculations (Ref. 12610 with a slightly different mixing param- 
eter, tan 2 #12 = 0.34 instead of 0.42 in the new calculations. 
However, note that the difference due to the parameter choice is 
quite small, as shown in the next subsection. 
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they included all the relevant neutrino processes in their 
calculations. 

Using the flux of each flavor neutrino on the Earth 
and cross sections of the relevant neutrino interactions at 
SK, we can calculate the expected event numbers from 
future galactic supernova neutrino bursts. SK is a water 
Cerenkov detector with 32 ktons of pure water, based at 
Kamioka in Japan. The relevant interactions of neutrinos 
with water are 



v e +p - 


-> e+- 


f n (CC), 


(11) 


V e + e~ - 


-> v e ^ 


-e~ (CC and NC), 


(12) 


v e + e~ - 


"> Ve - 


-e~ (CC and NC), 


(13) 


v x + eT - 


-> Vx " 


he" (NC), 


(14) 


V e +0 - 


-> e~ - 


f F (CC), 


(15) 


Ve+O - 


-> e+- 


f N (CC), 


(16) 



where CC and NC stand for charged- and neutral-current 
interactions, respectively. SK has restarted observation 
with lower performance, finishing the repair of the unfor- 
tunate accident which occurred on 12th November 2001. 
The effect of the accident on its performance is expected 
not to be serious for supernova neutrinos, because the 
fiducial volume does not change, and the threshold en- 
ergy change (from 5 MeV to about 7-8 MeV) influences 
the event number very little. Although the energy res- 
olution will become about \pl times worse, it does not 
matter for our considerations. In the calculations, we 
used the energy threshold and the energy resolution af- 
ter the accident, or in the so-called SK-II phase. For the 
cross section of the most relevant interaction Eq. (|llfl . 
we referred to the recent result by Strumia and Vissani 
|35j . and for the other processes we adopted the same 
ones as in Ref. [2(|. As a result, the expected events 
decreased by ~ 10% in contrast with the previous calcu- 
lation in Ref. [^fj] under the same conditions but for the 
v e p reaction cross section. 

The expected event number per unit energy range is 
shown in Figs. Eland[7| Figure [fj] shows the energy spec- 
trum of electrons (positrons) for the first 0.5 s with the 
LL20 model, which was obtained by conversion calcu- 
lation with the W02S and W95S presupernova models. 
When the magnetic field is not strong, B = 10 s G, the 
RSF conversions are absent, leading to an energy spec- 
trum independent of the presupernova model. (The slight 
difference between the two models comes from the differ- 
ence in the adopted mixing angle, i.e., tan 2 6*12 = 0.42 
for W02S and 0.34 for W95S.) On the other hand, when 
the magnetic field is sufficiently strong, Bq — 10 10 G, the 
energy spectrum is very sensitive to the adopted presu- 
pernova model, reflecting the results of conversion prob- 
abilities, Fig. \5\ Since the flavor conversions are not as 
efficient in the case of the calculations with the W02S 
model as those with the W95S model, the expected en- 
ergy spectrum with the W02S model is not as hard as 
that with the W95S model. 

Figure \7\ shows the energy spectrum at SK, obtained 
using the W02S and W02Z models with B = 10 8 , 10 10 




E e [MeV] 



FIG. 6: Energy spectrum of electrons (positrons) at SK for 
the first 0.5 s, obtained using the W02S and W95S presuper- 
nova models with Bo = 10 s , 10 10 G. The original spectrum of 
LL20 is adopted. 



G. Figure[7{a) is that for the LL20 model and Fig. Efb) is 
for the TBP20 model; the shapes of the spectra from the 
other TBP11 and TBP15 models are almost the same as 
that with the TBP20 model, although the absolute value 
is different. When Bq = 10 8 G, because the RSF conver- 
sions are absent and pure MSW flavor transitions occur, 
the energy spectra of the two models W02S and W02Z 
are degenerate. For the W02Z model with Bq = 10 10 
G, since efficient flavor conversions occur at low energies 
but not in the high energy region, the event numbers are 
suppressed only in the low energy region and high en- 
ergy tail is almost the same as that with Bo — 10 8 G. 
Since the average energy difference between the P e 's and 
j^.t's is not prominent for the TBP models, the differ- 
ence between energy spectra with various presupernova 
models is suppressed, compared with the LL20 model, 
particularly in the high energy tail. 



V. DISCUSSION 

A. What can we learn about the RSF effect from 
the neutrino signal? 

In the previous section, it was shown that the RSF 
conversion strongly depends on the presupernova models 
of various metallicities. The behavior of the flavor mix- 
ing in the supernova envelope is very different from one 
presupernova model to another, and they indicate various 
profiles in supernovae such as magnetic field strength and 
the value of (1 — 2Y e ) at the resonance points as well as 
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K e [MeV] 

FIG. 7: Energy spectrum of electrons (positrons) at SK ob- 
tained using the W02S and W02Z presupernova models with 
B = 10 s , 10 10 G. (a) The LL20 model for the first 0.5 s. (b) 
The TBP20 model for the first 0.25 s. 



the neutrino magnetic moment. Unfortunately, however, 
what we can observe is the energy spectrum of electrons 
(positrons) at detectors on the Earth alone, and much 
information on the detailed character is lost, and there- 
fore we can obtain only rather rough characteristics in 
principle. Here, we use as a simple indicator of the RSF 
conversions the following quantity: 

number of events for E e > 25 MeV . . 
number of events for E e < 20 MeV ' 

The values of Rsk for various models are summarized in 
Table [I] Even if we use the data for the first 0.25 or 0.5 
s, it appears that statistically sufficient discussions are 
possible. 

TABLE I: The values of Rsk for various models. Attached 
errors are statistical ones at the la level. 



B„ = 10 s G B = 10 10 G 

Model W02S W02Z W02S W02Z 

LL20 (0.5 s) 0.59 ± 0.02 0.59 ± 0.02 0.96 ± 0.04 0.72 ± 0.03 
LL20 (0.25 s) 0.41 ± 0.03 0.41 ± 0.03 0.71 ± 0.04 0.52 ± 0.03 
TBP20 0.25 ± 0.02 0.25 ± 0.02 0.35 ± 0.03 0.32 ± 0.03 

TBP15 0.21 ± 0.03 0.21 ± 0.03 0.33 ± 0.04 0.27 ± 0.03 

TBP11 0.22 ± 0.03 0.22 ± 0.03 0.35 ± 0.03 0.28 ± 0.03 



In practice, to make matters worse, there is a large un- 
certainty concerning the original neutrino spectrum emit- 
ted by the core collapse; actually the values of Rsk are 



very different between the LL20 and TBP20 models (see 
Table 0. Thus, we must reduce the systematic errors of 
the models in discussing the RSF effect from the energy 
spectrum obtained. At present, however, this problem 
is very difficult and there is no way but to wait for the 
future development of numerical simulation of supernova 
explosions. Therefore, at present, it is very difficult to say 
even whether the RSF effect actually occurred or not. 

If the systematic errors concerning the original neu- 
trino spectrum are considerably reduced by future de- 
velopment of numerical simulations, it is expected that 
useful implications for the RSF mechanism will be ob- 
tained from the value of Rsk- Therefore, at the end of 
this subsection, we discuss to what extent we can learn 
from the observed neutrino signal assuming that the sys- 
tematic errors are much reduced. If the magnetic field or 
the neutrino magnetic moment are too small to induce 
adiabatic RSF conversions, the metallicity of presuper- 
nova stars is not detectable as shown in the second and 
third columns of Table |l] On the other hand, in the case 
of a strong magnetic field like Bo = 10 10 G, the value 
of Rsk becomes larger than that for Bo = 10 8 G at the 
> (2-3)<7 level, which is a statistically meaningful num- 
ber. However, it is very difficult to estimate the metal- 
licity of the presupernova star, because the value of Rsk 
sensitively depends on the value of Bo, when it is larger 
than 10 9 G as shown in Ref. |2f|. Consequently, even 
if the systematic errors were considerably reduced in the 
future, all we could say concerning the RSF effect from 
the galactic supernova neutrino burst would be at most 
that the RSF conversions have occurred or not. Other 
detailed discussion concerning various quantities such as 
the magnetic field strength, supernova metallicity, and 
neutrino magnetic moment would not be accessible. 



B. Future prospects 

In the near future, it is expected that the many un- 
certainties which complicate the observations of the RSF 
conversion will be much reduced. These uncertainties 
are, e.g., the neutrino magnetic moment, the supernova 
magnetic field structure, and the metallicity of the su- 
pernova. 

If KamLAND receives a positive signal for the appear- 
ance of solar £> e , then it indicates that the spin- flavor 
conversion of the Majorana neutrino actually occurs in- 
side the Sun, and from the data we can obtain impli- 
cations for the nonzero value of the neutrino magnetic 
moment. Actually, in the Sun the relevant RSF conver- 
sion is v e <-> ^,,1-, whereas in supernovae it is v e <-> v^. T - 
However, if the neutrinos are Majorana particles, the ab- 
solute values of the transition magnetic moments, which 
are responsible for the above two processes, are the same 
(the sign is different; /iy = —[tji). Thus, depending on 
the value that is observed by the solar neutrinos, it can 
easily be estimated whether the RSF effect is the relevant 
conversion process for the supernova neutrinos; further, 
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if it is the relevant process, the magnetic field strength 
itself, not the combination /^„S, at the resonance points 
can be inferred. From future KamLAND results we can 
constrain only one of the magnetic moment tensor el- 
ements, which consist of three independent quantities, 
but we need all three values for the supernova case. Yet 
that information, if actually obtained, would be helpful 
for the estimation of other values and very useful. 

As for the metallicity, we will obtain useful inferences 
from astrophysical discussions. The lifetime of massive 
stars which end their life by gravitational collapse is much 
shorter than that of the Sun, and the progenitors of ob- 
served supernovae are, therefore, younger. Consequently, 
the metallicity of the galactic supernovae is expected to 
be at least the solar abundance or more metal rich. If the 
metallicity is higher than that of the Sun, the suppres- 
sion of (1 — 2 Yg) will be weaker and the RSF conversion 
will incline to be more adiabatic. On the other hand, the 
large and small Ma gella nic clouds are known to be very 
metal-poor systems |36|,|33. Thus, the determination of 
the explosion site will provide several rough estimations 
of the metallicity. Even though the explosion occurred 
in an optically thick environment such as the galactic 
center, the detected neutrinos alone could locate the su- 
pernova precisely; from the event number the distance to 
a supernova will be inferred, and from the event distri- 
bution the supernova direction can be determined |3S| . 

Another object that is related to the metallicity effect 
is the supernova relic neutrino (SRN) [39l Eof . Because 
the SRN is the accumulation of neutrinos from all the 
past supernovae, the SRN includes neutrinos from su- 
pernovae with quite low metallicity in the early phase 
of galaxy formation. Recently, the SK collaboration re- 
leased the first result of their search for the SRN sig- 
nal 0; an d h is only about a factor of 3 larger than 
the theoretical predictions by Ando et al. [39|. Data ac- 
cumulation for several more years may be decisive for 
actual detection and setting a severe constraint on the 
SRN flux. However, because the totally time-integrated 
neutrino spectrum during the neutrino burst is needed 
for the SRN calculation, it is very difficulty to estimate 
the SRN flux including the RSF effect, which is affected 
by the time-dependent density and magnetic field pro- 
files. Thus, it will be difficult to derive useful informa- 
tion about the RSF conversions from future SRN obser- 
vations. 



VI. CONCLUSION 

The RSF conversion of supernova neutrinos was inves- 
tigated in detail; this conversion is induced by a nonzero 



magnetic moment of neutrinos and a strong supernova 
magnetic field. Because of the effective potential for the 
RSF conversion which delicately depends on the devia- 
tion of Y e from 0.5, we studied the dependence of the RSF 
conversion on the adopted presupernova models. As the 
models, we used recent models by Woosley et al. [27j with 
solar and zero metallicities (W02S and W02Z) as well as 
that by Woosley and Weaver |2^| with solar metallicity 
(W95S), which was the only model adopted in our pre- 
vious paper 26], for comparison. 

The W02 models were calculated using a recent shell 
model and resulted in substantial revisions to the older 
data sets in the W95S model. As a result, for both 
W02 models the value of (1 — 2Y e ) suddenly drops at 
the radius where the RSF conversions occur, leading to 
a less efficient RSF conversion than that with the W95S 
model, whose (1 — 2Y e ) gradually decreases at the RSF 
regions. In particular, for the W02Z model, the value 
of (1 — 2Ye) becomes exactly zero at the boundary be- 
tween the Si + O and O + Ne + Mg layers, which makes 
the conversion almost completely nonadiabatic. Thus, 
we found that energy spectra expected at SK for the 
W02 models would not be as hard as those for the W95S 
model. Unfortunately, however, there are many uncer- 
tainties at present, which are concerned with, e.g., the 
original neutrino spectrum, the neutrino magnetic mo- 
ment, the supernova magnetic field structure, and the 
supernova metallicity. All these uncertainties complicate 
the discussions, and it is very difficult to obtain some use- 
ful information on the supernova metallicity or magnetic 
field strength, because the detected energy spectrum does 
not indicate which effect contributes to what extent. 

In the future, however, it is expected that the devel- 
opment of numerical simulations of supernova explosions 
will much reduce the uncertainties for the original neu- 
trino spectrum, and the KamLAND solar P e observation, 
whether detected or not, will provide useful information 
on the neutrino magnetic moment. We believe that all 
these developments will enable a more decisive investiga- 
tion of the RSF conversion of supernova neutrinos. 
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